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Universality of young cluster sequences 
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ABSTRACT 

Aims. Most stars do not form in isolation but as part of a cluster comprising anywhere between a few dozen to several 
million stars with stellar densities ranging from 0.01 to several 10^ MqPC~^. The majority of these clusters dissolve 
within 20 Myr. The general assumption is that clusters are born more or less over this entire density range. 
Methods. A new analysis of cluster observations is presented. 

Results. It demonstrates that, in fact, clustered star formation works under surprisingly tight constraints with respect 
to cluster size and density. 

Conclusions. The observed multitude of cluster densities simply results from snapshots of two sequences evolving in time 
along pre-defined tracks in the density-radius plane. This implies that the cluster size can actually be used to determine 
its age. 
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1. Introduction 



The very early phase of star formation is still poorly un- 
derstood because most stars form in dense clusters cm- 
bedded withi n giant molecular clouds (Lada & Lada 200l; 
iPudritz 2002[ ) often only observable at infrared wave- 
lengths. If the clusters consist of a thousand stars 
or more, they also contain massive O stars that 
drive the expulsion of gas from the cluster via stel- 
lar winds, ionisation and supernovae of early-type stars. 
This gas expulsion brings the clusters out of dynam- 
ical equilibrium, eventually exposing them and lead- 
ing to a cluster e xpansion w h ere the majority of stars 
become unbound ([H ills 1980; Goodwin fc Bastia n 2006t 
iBaumgardt fc Kroupa 2007; Bastian ct al. 2008), eventu- 
ally turning into field stars. Observations support this pic- 
ture because young 1 Myr) clusters are usually smaller 
than older (~ 20 Myr) clusters, but a systematic and quan- 
titative picture of the cluster expansion itself has so far been 
lacking. It is this evolution of exposed clusters over the first 
20 Myr that we address here. It is a time-span of special 
interest because it is precisely the period when planetary 
systems are developing. 

The observed densities of young star clusters range from 
less than 0.01 to several 10^ Mopc~^ as Fig. 1 shows (a list 
of cluster properties is given in Table 1 in the online ma- 
terial: Figer 2008, Wolff et al. 2007, Borissova et al. 2008), 
where the determination methods of the masses and radii 
from the literature sources are also explained). Here only 
clusters more massive than 1000 M© (massive clusters) are 
included, since the size and density of clusters containing 
only a few dozen stars are poorly determined. The wide 
variety of cluster densities has lead to the general assump- 
tion that clusters are also formed over this entire density 
range. As a result of the gas expulsion, a rapid dissolution 
of the cluster becaus e of their s uper- virial velocity disper- 
sion (H ills 1980; Kr oupa 2005t iGoodwin fc Bastian 20061: 
iGieles &: Bastian 20081 ) is then expected. 
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Fig. 1. Cluster density as a function of age for clusters more 
massive than 10"^ M©. The values were taken from Figer 
(2008), Wolff (2007), and Borissova et al. (2008) and refer- 
ences therein. See Table 1 in online material as well for a 
discussion of the error bars. 



2. Results 

When looking closely at the same clusters as in Fig. 1 in a 
different way by plotting instead the density as a function 
of the cluster radius R (definitions for the individual data 
points are given in the online material) leads to the remark- 
ably simple structure in Fig. 2. Apart from clusters existing 
in two groups, what this plot unambiguously shows is that 
these evolve along two well-defined tracks in the density- 
radius plane, strongly suggesting a bi-modal cluster evolu- 
tion. Moreover, as can be seen in Fig. 2, the two classes of 
clusters - in the following called starburst and leaky clus- 
ters (the terminology will be explained later) each start 
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from unique points in the mass-radius plane and develop 
at approximately the same speed along these two tracks. 
All clusters shown originate from one of these birth points 
in mass and radius. 

This is not so surprising for starburst clusters since their 
relatively similar mass has already led to speculations that 
they might have a common origin. For leaky clusters, how- 
ever, the prevailing view is that young clusters form more or 
less at arbitrary density and somehow dissolve afterwards 
on timescales of 20 Myrs. However, the new representa- 
tion of already known cluster parameters in Fig. 2 shows 
the common history of these clusters for the first time. The 
single track and sequential nature of the data points means 
that, contrary to the usual viewpoint, these clusters also 
start out with the same mass density and radial extent 
(and therefore mass) and all expand in the same way. 

Here a word about determinating of radii and densities 
in the different sources is nessesary. All starburst cluster 
data, apart from that for DBS 2003, were taken from Figer 
(2008). The data for DBS 2003 were taken from Borissova 
et al. (2008). In Figer (2008) the radius gives the average 
projected separation from the centroid position. The total 
cluster mass is given by extrapolating down to a lower mass 
cut-off of 1 solar mass, assuming a Salpeter initial mass 
function and an upper mass cut-off of 120 solar masses (ex- 
ceptions from this method are given in Figer 2008). The to- 
tal mass of the cluster was determined by assuming that the 
stars in the mass range 6-12 Mq constitute 5.5% of the total 
mass of the complete population of stars spanning the range 
0.1-100 M0. They compared their values with other mass- 
determining methods by Slesnick et al. (2002) and Weidner 
& Kroupa (2006) and found good agreement. In Borissova 
et al. the radius was defined by distance where the density 
profile exceeds twice the standard deviation of the surface 
density in the surrounding field. All values for the leaky 
clusters were taken from Wolff et al. (2007). There the ra- 
dius was determined as the median distance in degrees of 
the B stars in the sample from the centre of the cluster, so 
that, although not all cluster sizes were determined exactly 
in the same way, for each subset - the starburst clusters 
and the leaky clusters - the cluster radii were determined 
fairly consistently. The density was in all cases defined as 
p = 3Mc/4:TtR^, so it is extremely unlikely that the result 
is an artefact of the applied method. 

Quantitative analysis of Fig. 2 shows that, as the leaky 
clusters expand, they do not simply diffuse (which would 
give an i?~^-dependence), but the cluster density decreases 
as This means that the cluster loses mass via various 
processes such as outgassing, stellar mass loss, tidal evapo- 
ration, and escapers created through three-body encounters 
close to its centre. The relative importance of these mass- 
loss processes will need clarification in the future. As a di- 
rect consequence of the common history of these clusters, 
one can now determine the cluster size as a function of clus- 
ter age and finds a ''-dependence. As demonstrated 
in Fig. 3, the expansion velocity is ~ 2pc/Myr in the early 
phases and drops to ^Ipc/Myr at cluster ages '^20 Myr. 
In Fig. 3 it is shown that comparing simulations of clusters 
expandin g due to gas expulsion showing similar expansion 
velocity ( Olczak 20081 ) indicates that ~ 75% of the gas is 
expelled, though more detailed studies are required to con- 
firm this. However, this would already suggest that all leaky 
clusters shown in Fig. 2 formed with the approximately the 
same star formation efficiency of ^ 25%. 
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Fig. 2. Cluster density as a function of cluster size for clus- 
ters more massive than 10'^ Mq. The values were taken from 
Figer (2008), Wolff et al. (2007) and Borissova et al. (2008) 
and references therein. See although Table 1 in online ma- 
terial. 



It is thought that clusters and associations are formed 
with mass functions that are powerlaws with index -2. 
However, the results described here imply that the leaky 
clusters form in a narrow mass range (3.6 < log < 4.4) 
and develop quickly within 20 Myrs to lower masses (3.3 
< log Mc < 3.6). The small sample size currently does not 
allow further constriction of the primordial mass function 
for loss clusters. 

The situation is simpler for the starburst clusters. They 
all appear to be born with a cluster density of at least 10*^ 
M0pc~^ (red data points in the top left corner of Fig. 2) 
and then simply diffuse without further mass loss, i.e. a 
i?~'^dependence. In fact the expansion proceeds more or 
less linearly in time, i.e. as i? ~ ic (see Fig. 3). Having 
established a common history, one can easily deduce an ex- 
pansion velocity of 0.1-0.2 pc/Myr. From this expansion 
velocity, it should be possible to determine the star forma- 
tion efficiency through numerical simulations in the future. 
A possible reason for the expansion velocity in star-burst 
clusters to be lower than in leaky clusters is that the poten- 
tial well is too deep fo r the ionised gas to leave the cluster 
quickly (|Krouia 200l) . This means a much higher star for- 
mation efficiency for starburst clusters than for the massive 
clusters. 

Fig. 2 shows an additional feature, namely the location 
of the Galactic center clustre. Since the Galactic centre har- 
bours a massive black hole, one would expect its dynamics 
to differ considerably from starburst clusters. However, it 
lies directly on the starburst cluster track close to the data 
point of Arches. Given that the Galactic centre cluster is 
older than the Arches cluster, it follows that it must have 
started out at a higher initial mass and/or expanded more 
slowly. 

If the above systematics is true, then no clusters should 
exist with i? > 2 pc ages ages < 4Myr that possess masses 
higher than several 10'' Mq. An interesting consequence of 
the sequential nature of cluster development is that the ra- 
dius of the cluster can actually be used to determine its 
age. Cluster age determination is notoriously problematic. 
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Fig. 3. Cluster size as a function of cluster age for the same 
clusters in Fig. [21 The squares correspond to the leaky clus- 
ters and the circles to the starburst clusters. The dotted line 
indicates simulation results by Olczak (2009) for an ONC- 
like cluster with 75% gas expulsion. The drawn and the 
dashed lines are not fits but just there to guide the eye. 



striction that at least 200 cluster stars were detected. The 
embedded clusters are not located on the main evolution- 
ary sequence but form their own side-arm off the leaky clus- 
ter sequence below its extension to younger ages. It joins 
up when the embedded clusters become exposed after hav- 
ing reached a size of ~ 2-4 pc at an age of > 1 Myr. A 
further argument for the sequential nature of the cluster 
development is tha t the ONC is likely to be expanding 
()Kroupa et al. 200ll ). hence evolving along a time-density 
track. In other words, a relatively sparse cluster like Upper 
Cen Lup must once have been in a state similar to the 
Orion Nebular Cluster now the densest nearby cluster (see 
Fig. 4). Note as well that Ori lb, Ori Ic, and Ori la form a 
temporal sequence of once leaky clusters with initial densi- 
ties as high as the ONC. The initial stellar mass content of 
leaky clusters at the end of the embedded phase is at least 
3500 Mq. with a radius exceeding ~ 4pc, 

One would expect a similar side-arm for the starburst 
clusters, but so far no such embedded starburst clusters 
have been observed. However, this is not surprising since 
they would be very small {R < 0.1 pc) and located either 
close to the Galactic centre (at a distance of ~ 8kpc) or 
inside spiral arms within the Galactic plane. In addition, 
starburst clusters are intrinsically rare, and the embedded 
phase very short, so not many are expected to be found. 
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Fig. 4. Same as Fig. [2] but here the cluster density of em- 
bedded clusters with more than observed 200 members are 
also shown. The values were taken from Lada & Lada 2003. 

especially for very young clusters, where the uncertainty of 
the tracks of pre-main sequence star development can lead 
to very different age estimates. The density-radius relation 
defined here might be an easier and potentially more accu- 
rate method for cluster age determination. Clusters with a 
radius between 2-4 pc pose a special problem. Only the dis- 
tinction between very young leaky clusters (< 4Myr) and 
old starburst clusters (> 10 Myr) will require additional 
information due to their similar radii. However, as one can 
expect that in older starburst clusters some red supergiants 
are present, telling them apart is probably not too difficult. 

How do embedded clusters fit into this picture? Fig. 4 
again shows the cluster density as a function of size but also 
includes embedded clusters. The embedded clusters were 
selected from the data by Lada & Lada (2003) with the re- 



3. Discussion and conclusions 

The above results imply that star formation occurs only 
under an extremely limited set of conditions, and may 
therefore require a fundamental revision of its theories in 
our Galaxy, possibly including our own solar system. The 
"leaky clusters" could also be known as "mass-loss clus- 
ters", as this is what distinguishes them from starburst 
clusters." . For this leaky clus ter density regime, recent sim- 
ulations (jKiipper et al. 20"08l ) have found that the dynamics 
of embedded clusters show a strong correlation between the 
kinetic energy of the cluster stars and that of stars leaving 
the cluster, leading to a common sequence. This prediction 
might well be linked to the results described herein. 

The embedded cluster side-arm implies that many of 
the older clusters (like Ori la-c. Upper Sco, I Lac, etc.) 
must have gone through a phase like the ONC is in to- 
day. Thus the ONC can serve as model cluster for both 
the earlier phase of older leaky clusters and the future of 
the younger embedded clusters. Results for the ONC there- 
fore serve as a template for a multitude of other clusters. 
For example, observations and theory indicate that the cen- 
tral density in the ONC is so high that encounters between 
young disc-surro unded stars m ight affect t he protoplane - 
tary discs (Hcrbi g fc Terndrup 1986; Pfalz ner et al. 20"0^ 
and changes in the disc properties are likely to have direct 
consequences on the ongoing planet formation process. 

Finally the proposed reordering immediately raises the 
question of the origin of these two distinct cluster group- 
ings. The critical parameter is possibly the higher density, 
but could equally be the density-temperature ratio or dif- 
ferent types of turbulence, etc. The obvious task for the 
near future is to understand why these specific conditions 
are favoured in our Galaxy and whether the same applies 
in other spiral galaxies, too. 
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Table 1. P roperties of clusters an d a ssociations more m assive than 10"^ in the Galaxy ^ (jFiger 20081 and references 
therein), ^ (jBorissova et al. 20081 ). ^ (j Wolff et al. 200"7l and references therein). The errors in age result from (a) their 
deduction from theoretical pre-main-sequence star tracks for young stars, which show large discrepancies; (b) the age of 
the cluster is found by averaging over the so deduced ages of the single stars; and (c) the age determination depends on 
the knowledge of the distance to the cluster - a parameter that is often poorly constrained. In cases where there is no 
error bar this does not mean that the age is accurately known, but rather that the authors did not give any error value. 
Usually there are no error values quoted for the mass or the density of the cluster, nor is the cluster size consistently 
defined for the different clusters. An error of a factor of 2 is probably realistic and can be expected to be even larger for 
the more distant clusters, so for this and the following figures one should always keep in mind that the errors of these 
observed values are rather large. 
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T able 2. Properties o f embedded clusters containing more than 200 observed stars. 
^ ()Lada fc Lada 20^ . and references therein) 
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